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Hyperfine interaction and magnetoresistance in organic semiconductors 
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We explore the possibility that hyperfine interaction causes the recently discovered organic mag- 
netoresistance (OMAR) effect. Our study employs both experiment and theoretical modelling. An 
excitonic pair mechanism model based on hyperfine interaction, previously suggested by others to 
explain magnetic field effects in organics, is examined. Whereas this model can explain a few key 
aspects of the experimental data, we, however, uncover several fundamental contradictions as well. 
By varying the injection efficiency for minority carriers in the devices, we show experimentally that 
OMAR is only weakly dependent on the ratio between excitons formed and carriers injected, likely 
excluding any excitonic effect as the origin of OMAR. 

PACS numbers: 73.50.-h,73.50.Qt, 



I. INTRODUCTION 

Organic 7r-conjugated materials have been used to 
manufacture devices such as organic light-emitting diodes 
(OLEDs) [ll|L photovoltaic cells 0,110 and field-effect 
transistors [y, BIB - Recently there has been growing in- 
terest in spi n lid ITTl If2l and magnetic field effects 
(MFE) [ll llllllllJll 13 IM El JPJP in 
these materials. Frankevich and coworkers |24l |25| and 
Kalinowski and coworkers [22l [2(J studied the effect of 
small magnetic fields on excitonic processes that occur in 
OLEDs, such as photoconductivity, electroluminescence 
(EL) and exciton dissociation at the electrodes. They 
have convincingly explained their findings using a model 
where the applied magnetic field reduces the effect of the 
hyperfine interaction between electron/hole spin and the 
hydrogen nuclei in the organic molecules. We will refer 
to this work as the excitonic pair mechanism model that 
will be treated in detail in section If V CI 

We recently discovered 01 a large and intriguing mag- 
netoresistive effect in OLEDs, which we dubbed organic 
magnetoresistance (OMAR). OMAR may find applica- 
tion in magnetic field sensors, e.g. in OLED interactive 
displays (patent pending, see demonstration video in Ref. 
|la|). In addition to its potential applications, OMAR 
poses a significant scientific puzzle since it is, to the best 
of our knowledge, the only known example of large room 
temperature magnetoresistance in non-magnetic materi- 
als with the exception of narrow-gap high-mobility ma- 
terials [27| . 

In the present work we examine whether the excitonic 
pair mechanism model [2^, 0, Hi| can explain OMAR. 
For this purpose we recast this model (see section ITV Cf) 
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FIG. 1: A schematic drawing of the device and the magne- 
toresistance experiment 



into a form suitable for discussing the whole body of ex- 
perimental MFE data, including both EL and transport 
measurements (OMAR). We find that whereas the ex- 
citonic pair mechanism model can explain several key 
aspects of experiments, it leads to several serious con- 
tradictions with experiment, especially in relation to the 
magnetotransport data. We trace the origin of these con- 
tradictions to the fact that this model in its present form 
has been applied to describe the spin-dynamics of neutral 
polaron pairs that do not significantly affect the current. 
Furthermore, we show experimentally that the magni- 
tude of OMAR depends only weakly on the ratio, r/i , be- 
tween excitons formed and charge carriers injected into 
the devices. In our opinion, this shows that any model 
based on excitonic processes fails to explain OMAR. 



II. EXPERIMENTAL 

We first describe the sources for the various organic 
semiconductors we used in our study. The 7r-conjugated 
polymer polyfluorene (PFO) was purchased from Amer- 
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ican Dye Source, inc. The methyl-substituted ladder- 
type poly(p-phcnylcne) (McLPPP) polymer was synthe- 
sized as described elsewhere [2£j. The 7r-conjugated small 
molecule Alq3 was purchased from H. W. Sands corpo- 
ration and was used as received. 

The fabrication of the organic sandwich devices started 
with glass substrates coated with 40nm of indium-tin- 
oxide (ITO), purchased from Delta Technologies. The 
conducting polymer Poly (3,4-ethylenedioxythiophene)- 
poly (styrenesulfonate) (PEDOT), purchased from H. C. 
Starck was spin coated at 2000 rpm on top of the ITO 
to provide an efficient hole injecting electrode. All other 
manufacturing steps were carried out in a nitrogen glove 
box. The active polymer film was spin coated onto the 
substrate from a chloroform solution. The small molec- 
ular film layers were made by thermal evaporation. The 
cathode, either Ca (with an Al capping layer), Al, or Au, 
was then deposited by thermal (Ca) or electron beam 
evaporation (Al, Au) at a base pressure of « 1 x 10 -6 
mbar on top of the organic thin films. The device area 
was 1mm 2 for all devices. The general device structure 
used for our measurements was metal/organic semicon- 
ductor/metal (see Fig. Pi. 

A schematic drawing of the device and the experiment 
are shown in Fig. ^ The samples were mounted on the 
cold finger of a closed-cycle helium cryostat located be- 
tween the poles of an electromagnet. The magnetocon- 
ductance ratio, AI/I, was determined by measuring the 
current, /, at a constant applied voltage, V, for different 
magnetic fields, B. EL of the devices was measured with 
a photomultiplier tube that was shielded from the mag- 
netic field using high-saturation mu-shield foil. All the 
reported data are for room temperature. 



III. EXPERIMENTAL RESULTS, PART I 

A. Magnetic Field Effect on Current and 
Electroluminescence 

OMAR devices have the unique property of showing 
large magnetoresistance while being also highly electro- 
luminescent. Fig.EJand Fig.|3|show AI/I and AEL/EL 
vs. B in a PEDOT/PFO/Ca and PEDOT/Alq 3 /Ca de- 
vice, respectively, measured at a constant voltage. These 
data show that the MFE exists both in the electric and 
luminescent measurements with comparable magnitude. 
Note that the shape of (AI/I)(B) and (AEL/EL)(B) 
are equivalent and that both scale in the same manner 
upon changing V. Both effects therefore share a common 
origin. EL and / are related through: 

EL cx r/I, (1) 

where r\ is the EL quantum efficiency. For the sake of 
the following discussion it is useful to break down the 
EL process into three steps , and r\ can accordingly be 
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FIG. 2: Magnetic field effect (MFE) on current (bold) and EL 
in a PEDOT/PFO(« 100 nm)/Ca device measured at several 
different constant voltages at room temperature. The MFE 
on EL measured at a constant current, 5 [ih corresponding 
to ~ 13 V, is also shown. 
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FIG. 3: MFE on current (bold) and EL in a PEDOT/Alq 3 
(~100nm)/Ca device measured at several different constant 
voltages at room temperature. The MFE on EL measured at 
a constant current, 11 //A corresponding to « 10V^, is also 
shown. 



TABLE I: Function and parameter values used for fits in Fig. [I] 



fitting function: B' 2 /(B' A + Bg) 


material name 


Bo/mT 


regio-regular P3HT 


5.1 


pentacene 




5.8 


Weller et al. 1 




5.7 


Weller et al. 2 




13.7 


fitting function: (B/(|B| + B )f 


material name 


B /mT 


PFO 




5.4 


Alq 3 




5.4 


PPE 




5.5 


Yoshida et al. 




14 



B(mT) 

-33.3 0.0 33.3 66.7 100.0 133.3 166.7 200.0 
1 6 -I < I < 1 < 1 > 1 < 1 < 1 < 1 



written as 



i] = 771/72773, 



(2) 



where 771 is the fraction of the injected carriers that form 
electron-hole pairs, 772 is the fraction of the total number 
of excitons that are singlets, and 773 is the singlet emission 
quantum efficiency. 

There are two distinct possibilities regarding a poten- 
tial explanation of the MFE on EL (MEL): (i) I = I{B) 
and (ii) rj = 7/(73) . In terms of the underlying mechanism, 
a type (i) model involves the charge carriers directly, 
whereas a type (ii) model is based on excitonic mech- 
anisms. Frankevich and coworkers |24| and Kalinowski 
and coworkers (2^, have previously used the excitonic 
pair mechanism model (i.e. a type (ii) model) for ex- 
plaining their findings regarding the MEL effect. The 
question naturally arises whether their model can also 
explain OMAR. To answer this question we recast the 
pair-mechanism model in a form suitable to simultane- 
ously examine the prediction of this model for AEL/EL 
and Al /I. Before proceeding, however, we want to sum- 
marize the experimental evidence that supports the claim 
that the hyperfine interaction causes MEL and OMAR. 



IV. DISCUSSION, PART I 

Fig. 21 shows a summary of measured OMAR traces, 
taken from our previous publication 23], in PE- 
DOT/Organic Layer/Ca devices employing several dif- 
ferent organic layers. The traces of the various data 
sets have been normalized to achieve a suitable graph- 
ical representation. The solid curves are fits using em- 
pirical laws of the forms (B/(\B\ + Bq)) 2 (thicker lines) 
and B 2 /(B 2 + Bq) (thinner lines); the values for the fit- 
ting parameter, Bq, are given in table[I] These empirical 
laws were introduced by us in a previous publication |23J 
and yield excellent one-parameter fits. For comparison, 
data on the MEL effect in a poly(phenylene-vinylenene) 
OLED measured by Yoshida et al. [Tg an d that on 
the triplet photoyield in organic solutions measured by 
Weller et al. I29t l30l are also shown. Furthermore it is 
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FIG. 4: Normalized OMAR traces in PEDOT/Organic 
Layer/Ca devices, with PFO, Akj3, poly(phenylene- 
ethynelene) (PPE) , pentacene and regio-regular P3HT as the 
organic layers taken from Ref. |23H . For comparison, data on 
the MEL effect in a p oly ( pheny lene- viny lene) OLED mea- 
sured by Yoshida et al. flq| and that on the triplet pho- 
toyield in solutions of organic materials measured by Weller 
et al. |29l |30| are also shown. The data are shown as scatter 
plots as detailed in the legend. The solid curves are fits using 
empirical laws of the forms (B/(\B\ + Bo)) 2 (thicker lines) 
and B 2 /(B 2 + Bq) (thinner lines). Please note that the data 
by Yoshida et al. and the corresponding fit refer to the upper 
x-axis, whereas all other data refer to the lower x-axis. 



shown that the data by Yoshida et al. and Weller et al. 
can also be fitted accurately by our empirical laws sug- 
gesting a common origin of OMAR and other MFE data. 
Since Yoshida et al. and Weller et al. have interpreted 
their data in terms of the pair mechanism (hyperfine in- 
teraction) model, this suggests that OMAR is also caused 
by hyperfine interaction. As a matter of fact, Weller 
and coworkers [2^, |3(} have directly shown through tran- 
sient measurements of the delayed fluorescence caused by 
triplet-triplet annihilation in solutions of organic mate- 
rials that a weak magnetic field can modulate the spin 
multiplicity within 10-20ns after optical excitation. Fur- 
thermore, Schulten et al. jU showed that the width 
of the MFE traces can be calculated from first princi- 
ples based on the hyperfine coupling constants. However, 
Schulten et al. 31] did not specify an analytical result 
for the dependence of the MFE on B and the origin of the 
simple, analytical fitting formulas, (B/(\B\ + Bq)) 2 and 
B 2 /(B 2 + B 2 ), is not yet well-established. We have there- 
fore performed a simplified calculation of the dependence 
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of the MFE on B (see Appendix) where we show that the 
origin of the empirical law, B 2 /(B 2 + B 2 ), can be readily 
understood. We note that this formula is closely related 
to the Lorentzian function B 2 /(B 2 + B 2 ). However, we 
do not yet understand the origin of the other empirical 
law, B 2 /(\B\+B ) 2 . 

A. What determines the width of OMAR traces? 

The simple model calculation presented in the Ap- 

1/2 

pendix gives the relation Bq = (J2i a H %) f° r tne width 
of the OMAR traces, whereas Schulten's theory [3lJ gives 

1/2 

Bo = v3 (J2i a H i) f° r nuclear spin 1/2. In either 
case, it should therefore be possible to calculate Bg from 
published values for the hyperfine splitting, an in the 
electron spin resonance spectra of the organic molecules 
(see e.g. Ref. As a matter of fact, the width of 

the Weller et al. [2|j data numerically coincides with 
Schulten's formula with pretty good accuracy. However, 
for our OMAR data the experimental width is consid- 
erably greater than the calculated one. For pentacene, 
for example, we calculate Bq = lmT and Bq = 1.8mT 
using our and Schulten's formula, respectively. The ex- 
perimental value is however Bq = 5.8rnT. Indeed, the 
fact that the widths of the MFE traces in OLEDs are 
much wider than expected has already been recognized 
by Yoshida et al. 18]. Yoshida et al. therefore sug- 
gested that this is a result of lifetime broadening (see 
section llV Cjl . Our experiments however exclude this pos- 
sibility: The pair lifetime should be a sensitive function of 
temperature and should vary considerably between differ- 
ent materials in contradiction with experiment (see Fig. 01 
and Ref. [23j|). Although we have been unable to find a 
convincing explanation for this discrepancy, it is possible 
that hyperfine coupling with hydrogen nuclei of neighbor- 
ing molecules in the densely packed films and anisotropic 
hyperfine coupling may account for (part of) the excess 
width. A further disagreement between theoretical ex- 
pectation and experimental results is evident from the 
observation that Bq is similar in both small molecules 
and polymers. This is unexpected because McConnell's 
relationship |33| states that an.i — Qpi, where pt is the 
spin density at nucleus, i, and Q ps 3mT for conjugated 
molecules |32j . For a polymer with N repeat units we 
therefore have 

a H cx N~ 1/2 , (4) 

where ^ denotes a sum over all nuclei in the polymer, 
whereas is a sum only over a single repeat unit, pi is 
the spin-density at nucleus, i, whereas p\ is the spin den- 
sity at nucleus, i, in the corresponding monomer. This 
result implies that the MFE cones should be consider- 
ably narrower in polymers than in small molecules such 



as Alq3, in contradiction with the experimental results 
(see Fig.SJl. 

B. Potential implication of hyperfine coupling to 
organic spintronics application 

In the preceding paragraphs we have summarized some 
of the experimental evidence for the importance of hyper- 
fine coupling in organics. Furthermore, the discovery of 
OMAR and its possible relation to hyperfine interaction 
illustrates the importance of the study of spin-dynamics 
in relation to transport phenomena. In particular, we 
believe that the existence of hyperfine interaction has 
fundamental implications to the currently emerging field 
of organic spintronics |ld 1 1 1| . Since hyperfine coupling 
leads to time-evolution of the electron spin and since 
the local nuclear spin configuration is different for each 
molecule in the film, hyperfine coupling will lead to spin- 
decoherence. The relevant decoherence time-scale, T2 is 
given by 

T 2 ps 9 -^Bq ps Ins (5) 

Assuming the drift term to be dominant in OLEDs, we 
therefore obtain for the spin-transport length, A 

A ps uFT 2 ps lA (6) 

where we have used p = 10~ 4 cm 2 (Vs)' 1 for a typical 
mobility value and F = lQ^Vcm^ 1 for a typical value 
for the electrical field. This result clearly shows that hy- 
perfine interaction may seriously limit the spin-transport 
properties of organic semiconductors. Of course, the ef- 
fect of hyperfine coupling can be switched off by applying 
B > Bq, a property which may have some use in spin- 
tronics applications. 

C. Excitonic pair mechanism model 

Next, we examine whether OMAR can be explained 
by the excitonic pair mechanism model that was applied 
to MFE in OLEDs by Frankevich |24[ and Kalinowski 
|22^. When electrons and holes are injected from the 
cathode and anode into the organic laye r, they form neg- 
ative and positive polarons, respectively [34|, |35|]. As long 
as the distance between positive and negative carriers is 
larger than the Coulomb capture radius, rc, they do not 
feel each other's attraction, and we refer to them as free 
charges. In order for the pair mechanism model [22, 13] 
to work we assume that at separations less than rc, the 
carriers are organized in bound pairs, which we refer to as 
a polaron pair (PP). At separations less than the single 
particle wavefunction extend, the exchange interaction 
becomes important and the pair has to be represented 
by a single, properly symmetrized, wavefunction. We 
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FIG. 5: The schematic energy level diagram illustrates the 
simplest possible pair mechanism model which includes three 
different species: (i) free charges with population C, (ii) po- 
laron pairs PPs, and (iii) singlet Se and triplet excitons Te 
with large exchange interaction, J. The various transition 
rates are indicated. 



refer to this as an exciton. We note that the pair mecha- 
nism model is an example of a spin-dependent effect that 
does not require (thermal) spin-polarization. One way 
of understanding this is to realize that, since the carriers 
form pairs, the "ensemble" consists of two spins only and 
very large "effective spin-polarization" is automatically 
achieved. The formation of pairs is therefore essential to 
this model. 

Closely following the treatment by Frankevich [24[ we 
now formulate rate equations that govern the population 
of free charges, polaron pairs and singlet and triplet exci- 
tons. The relevant levels and transition rates are shown 
in Fig. 1 PP, 3 PP , 3 PP+, 3 PP- denote the pair 
populations, where the superscript denotes the multiplic- 
ity, whereas the subscript denotes the spin-projection. 
Because of negligible exchange interaction in PP, the 
four PP states are degenerate for zero applied B-field. 
The basic idea of the pair mechanism model is that the 
PPs may under go a spin motion induced by the hyper- 
fine interaction |3lt 13a . I37L l3q (see Appendix) . We de- 
note the rate of conversion between 1 PP and 3 PPo as 
kHF <n and that between other PP as Uhf- If B=0, then 
kuF = kuFa- An applied B-field leads to Zeeman split- 
ting, AE between levels of different spin-projection. If 
AE > HkuFoi then Jzhf = 0, however, kup remains 
unchanged. In summary of this paragraph, the simplest 
possible model for explaining the magnetic field effect of I 
and EL requires rate equations for three different species, 
namely (i) free charges that carry the current, (ii) pairs, 
such that it becomes meaningful to talk about singlet 
and triplet states, and (iii) excitons that, by virtue of 
their large exchange energy, furnish the model with spin- 
dependent pair recombination rates. We note that pairs 
and excitons do not contribute to the current, because of 
their overall neutrality. The various rate equations are 
given by: 



Gc + k P c(j^{ i PP))-kcpC = Q (7) 

-k CP C + k HFo ( 3 PP - 1 PP) + k HF { 3 PP+ - 1 PP) 

+k HF ( 3 PP- - 1 PP) - (kp C + kpsfPP = (8) 

X -k CP C + k HFo CPP - 3 PPo) + k HF ( 3 PP+ - 3 PP ) 

+k HF ( 3 PP_ - 3 PP ) - (k PC + k PT ) 3 PP a = (9) 

-k CP C + k HF CPP - 3 PP ± ) + k HF { 3 PP - 3 PP±) 

+k HF ( 3 PP T - 3 PP±) - (k PC + k PT ) 3 PP ± = (10) 

The definitions of the various transition rates are given 
in Fig. The first equation is the rate equation for 
the free charges, C. Gc is the generation rate for C, 
which is equal to the rate of carrier injection minus 
the rate of emission of carriers at the electrodes. The 
next three equations are the spin-dependent rate equa- 
tions for PPs. We neglected the upwards transitions 
ksp,kTP,kcs (dashed arrows in Fig. [SJ) to simplify the 
rate equations. We assume that Gc is constant as long 
as the voltage remains fixed, in particular it is indepen- 
dent of B [23( . We obtain the following solutions to the 
rate equations: 



AI 
AEL 



kpc 



= m 



(l-r) 2 



(4fe5+r + 3)[fe§(r + 3) + 2(r + i)] 



EL 



1 -r 



f^(r + 3) + 2(r + l) 



R = 



AI 
I 



AEL 



EL 



m- 



kpc 



AEL 



EL 



v 

1 - r 



2[2|^+r + l] 



(12) 

(13) 
(14) 



where we have used EL cx 1 PP and / cx C; r = 
kps/kpT 0. The notations |y and |/ refer to mea- 
surements where the voltage, V, or the current, I, were 
kept constant, respectively. We have used kHF ^> 
kps,kpT,kpc (see section llV A|l . 

We will now show that these results are in clear contra- 
diction with experiments. Equ. (|llf) shows that AI / I\v 
is always positive in this model in contradiction with the 
experimental results where both positive and negative 
magnetoconductance is observed |23J. Furthermore, it 
is seen that AI/I\y is a second order effect, whereas 
AEL I EL\y appears in first order in (1 — r). This con- 
tradicts the experimental observation (see Figs. and 0J) 
that there exists a fixed ratio, R, between the two effects. 
The model result for R (Equ. I|13fl ). however, is a func- 
tion of r and kpc/kpx- In the present case of polyflu- 
orene and Alq 3 OLEDs, the positive sign of R can be 
brought into agreement with the model only if r < 1. So 
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in the following discussion we will limit ourselves to this 
case. For the sake of discussion, we will study two limits: 
kpc/kpT > 1 and kp C /k PT < 1. If kp C /k PT > 1, Jill) 
(|12fl H14fl can be simplified as: 



AI 

T 



AEL 



4fe(r + 3)' EL 



kpc ( r 



AEL 



EL 



3) 
1 - r 

^ kpc 
kpT 



(15) 



If kpc/kp T < 1, (jnj lO JUJ can be simplified as: 



A/ 



kpc 
kpT 



{l-rf 



AEL 



2(r + l)(r + 3)' EL 
AEL 



EL 



1 - r 
2(r + 1) 

1- r 
2(r + 1) 



(16) 



From these results it follows that %^| v 



can be large 

only if kpc ~ &pt- Whereas this condition could be acci- 
dentally satisfied in a single material at a single temper- 
ature, it cannot hold true at all temperatures in different 
materials. The model's conclusion that =f-\ v is always 
small results because the free charges, the only current 
carrying species in the model, do not directly participate 
in the spin-dependent reactions or the hyperfine transi- 
tions. In particular, when kpc/kpT 3> 1 most polarons 
exist as free charges and the pair concentration and there- 
fore the MFE are small. When kpc/kpx <C 1 the pair 
concentration and therefore the MFE are large, but the 
large change in the pair population does not affect the 
current because the pairs do not dissociate. Equs. i|15|) 
and i|16|) highlight another contradiction between model 
and experiment: In the present model it makes little dif- 
ference whether AEL /EL is measured while the voltage 
or current are kept constant, whereas the experimental 
results (Figs. |5] and |3J showed that AEL / EL\v is about 
ten times bigger than AEL/EL\j. 



V. EXPERIMENTAL RESULTS, PART II 

In this section we determine experimentally whether 
OMAR is related to an excitonic effect or not. This is 
possible, since 771 that appears in Equ. 1|12|) can be var- 
ied experimentally by several orders of magnitude. By 
varying the corresponding electrode materials one can 
control the injection of minority charge carriers, so that 
the exciton formation can be tuned with ideally minimal 
effect on the current density, which is determined mostly 
by the majority carriers. This can be easily achieved in 
hole-dominated PFO devices by choosing cathode (top 
electrode) materials with different work functions, but 
it is difficult in practice to fabricate electron-only Alq3 
devices with different anode (bottom electrode) materi- 
als, because efficient electron injection requires reactive 
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FIG. 6: a) Magnitude of the magnetoconductance, AI//, in 
PFO (see inset) as a function of the exicton/carrier ratio 771 
at B = lOOmT. b) EL as a function of current, c) Current- 
voltage (I-V) characteristics. All data were obtained at room 
temperature. 



metals such as Ca. We found, e.g. in a Ca/Alq3/Ca de- 
vice, that the bottom Ca electrode oxidizes quickly before 
and/or after evaporating Alq3 on top of it. So, instead of 
Alq3 devices we studied a second polymer, MeLPPP (see 
Fig. inset) in addition to PFO, to show that the con- 
clusions we draw are not limited to a particular choice of 
polymer. MeLPPP is a suitable choice, because it shows, 
like PFO devices, large OMAR as well as intense EL. 

Three types each of PFO and MeLPPP devices were 
fabricated and measured using Ca, Al and Au as cath- 
ode materials. We note that whereas only one set of data 
points for each device is shown, the reported experiments 
were repeated several times and very reproducible results 
were obtained. Figs. El and [7| show the current- voltage (I- 
V) characteristics (Figs. anddO, the measured EL 
intensity as a function of I (Figs. EJ) and [3)), and the 
magnitude of AI /I as a function of rji (Figs. EJi, and|J]i). 
The exciton/carrier ratio, 771 was inferred from the data 
shown in panels b) where we found that the magnitude of 
EL, at a given current, in Ca devices is about one order 
of magnitude larger than that in Al devices, and nearly 
three orders of magnitude larger than in Au devices. This 
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lated to an excitonic effect, but caused by a direct effect 
of B on the (mobility of the) free carriers, we would have 
expected no dependence at all, which does not match the 
measurement either. However, the observed weak depen- 
dence can easily be attributed to several possible factors: 
It is possible that the interface resistance of polymer/Au 
is larger than that of polymer/Ca resulting in additional 
resistance that is not subject to OMAR. Previous stud- 
ies |40| report that Au indeed leads to a non-Ohmic top 
contact, possibly because the wetting of Au and there- 
fore the physical contact is inferior, or because Au de- 
position, which has to be evaporated at a much higher 
temperature than Ca, leads to damage of the underlying 
polymer surface. Moreover, in unipolar devices space- 
charge limited current conditions occur which are possi- 
bly unfavorable for OMAR. In bipolar devices, however, 
the space charge of the two carrier types partially cancel 
each other. The fact that changing the cathode material 
leads to additional effects, rather than merely changing 
771, is shown in Figs. [Ht and[7|; where it is seen that Au 
cathode devices show a significantly increased device re- 
sistance, probably due to the above mentioned parasitic 
resistance contributions. We have however been unable 
to come up with any probable or improbable explanation 
how the dependence on 771 could be sub-linear in an ex- 
citonic model. We could imagine superlinear behavior, if 
the spin-dependent reactions are bimolecular. Our mea- 
surements therefore show that OMAR is most likely not 
related to an excitonic effect. 



FIG. 7: a) Magnitude of the magnetoconductance, AI/I, in 
MeLPPP (see inset) as a function of the exicton/carrier ratio 
r\\ at B — lOOmT. b) EL as a function of current, c) Current- 
voltage (I-V) characteristics. All data were obtained at room 
temperature. 



is evidence of the increasing difficulty of electron injec- 
tion from the cathode into the polymer going from low 
work-function Ca to high work-function Au. This is well- 
known 0] to result from the mismatch of the cathode 
work function and the polymer's conduction band. Cor- 
respondingly, rj is one (three) orders of magnitude lower 
in Al (Au) devices compared to Ca devices. Since, 772 and 
773 in Equ. @ are properties of the excitons [39( in each 
organic material, and they do therefore not depend on 
the carrier injection efficiency (hence the cathode mate- 
rials), we simply use 77 as a measure of the exciton/carrier 
ratio, ?/i , of the various devices. 

Most importantly, Figs.Et an d show that the mag- 
nitude of AI/I increases as 771 increases in both PFO and 
MeLPPP devices (results for three choices of the current 
are shown for each device) . At first sight this trend seems 
to confirm an excitonic origin of OMAR. However, closer 
inspection shows that the 771 dependence is much weaker 
than expected. Whereas a linear dependence on 771 is ex- 
pected for an excitonic effect, we find AI/I cx 77", with a 
ranging from 1/3 to 1/2. If, however, OMAR is not re- 



VI. DISCUSSION, PART II 

The excitonic pair mechanism model fails to explain 
the comprehensive experimental data, but it still suc- 
ceeds to account for several experimental observations 
(namely the magnitude and "universality" of Bq) that 
would otherwise be hard to explain with any model that 
does not involve the nuclear spin. So, it appears that a 
model is demanded, which retains some of the features of 
this model, namely the hyperfine interaction coupled to 
spin-dependent interactions, but applies them to a spin- 
dependent reaction that directly influences the mobil- 
ity. The suggestion by Frankevi ch 12m to consider spin- 
dependent bipolaron formation [34j appears promising. 
In the following we will demonstrate that the assumption 
that B acts directly on the mobility, /1, correctly repro- 
duces the experimental data. We are however presently 
not able to identify or formulate a specific mechanism 
with the required properties. Such a theory lies far be- 
yond the scope of the present paper and we hope that 
the scientific community will attempt to formulate such 
a theory. 

We now wish to discuss the dependence of EL and I 
on the mobility. The most general formula for for the 
current density, </, is: 

J = e(p/j p - 77/i„)F (17) 
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where e is the elementary charge, p and n are the den- 
sity of positive and negative polarons, respectively and 
\x p and [i n are their respective mobilities, F is the elec- 
tric field. All these quantities may in general depend 
on the position, x. In the case of unipolar transport, J 
in OLEDs is often assumed to be given by in the Space 
Charge Limited Current (SCLC) formula: 



In summary, a pair mechanism model that acts di- 
rectly on the current carrying species naturally implies 
the observed intimate relationship between AEL/ EL 
and AI /I. Furthermore it implies the experimental ob- 
servation that AEL/EL\i < AEL/EL\ V (see Figs. 
andEJl. 



9 V 



(18) 



where V is the applied voltage and d is the device thick- 
ness, e and eo have their usual meaning. This formula 
holds true only in the trap-free limit, in the case of traps 
the term CK) in Equ. (|18|1 must be replaced by a more 
general function, / [41|. In the following we will 
assume that the mobility of one type of carrier greatly 
exceeds that of the other, specifically fi p 3> fi n - This 
assumption is usually fulfilled in OLEDs. 

The relation for EL depends on the the relative mag- 
nitude of recombination time (r rec ) and transit time {jt) 
|42|. The ratio r rec / (r rec + r t ) is a measure of the proba- 
bility that the carrier will transit the device without re- 
combining with another carrier. For injection-controlled 
EL (jt < r rec ) we obtain: 



EL oc 7 



Jtotal 

AEL 
EL 



Jnjp 

Op 

A7 _ 
T ~ 



Hp 



(19) 
(20) 
(21) 



7 is the electron- hole recombination constant, which is 
usually assumed to be given by the Langevin recombina- 
tion formula 0, 7 = e(p n + fi p )/ee . Equ. (|2*T|) shows 
that the experimentally observed intimate relation be- 



tween 



AEL 



, and ^M- naturally follows in a type (i) model. 
However, result Equ. (|21|l implies that the two effects 
should be strictly of equal magnitude, which is not quite 
the case experimentally (see Figs. El and |3J) where AJLL 
is up to twice as big as . A more general relationship 
is therefore required. For the case r™ > r rec > rf we 
obtain: 



EL cx j n 



Jtotal 

AEL 
EL 



Jp 

A^ r , 
Hn 



AJ 
T 



Hp 



Furthermore, for volume controlled EL (r t > 77 
/!„, we obtain: 



which requires fi p 

AEL 
EL 

AI 
T 



Afj, p 



+ 



(m p ) 



A/ip 



Afi„ 
Afi p 



Afi p 



fJ.p 



(22) 
(23) 

(24) 
0), 



(25) 
(26) 



VII. CONCLUSION 

We have explored the possibility that hyperfine inter- 
action causes the recently discovered organic magnetore- 
sistance effect using both experimental data and theo- 
retical models. Other groups have previously explained 
magnetic field effects on excitonic processes in organics 
using the hyperfine pair mechanism model. We show that 
both kinds of data can be fitted using the same empirical 
laws, either of the form B 2 /(B 2 + B%) or B 2 /(\B\ + B ) 2 
dependent on material. The only fitting parameter, Bq, 
assumes values that, at least at first sight, seem typical 
of hyperfine interaction. This suggests that OMAR is 
caused by hyperfine interaction as well. We succeeded to 
"derive" the empirical law B 2 /(B 2 + Bq) from the stan- 
dard hyperfine Hamiltonian. We also showed that hy- 
perfine interaction may seriously limit the spin transport 
length, which is of primary importance in spintronics ap- 
plications. 

In order to further test the hyperfine interaction hy- 
pothesis, we examined a pair mechanism model, sug- 
gested previously by other authors which considers the 
electron spin precession caused by the hydrogen nuclear 
magnetic field. However, we found several fundamental 
contradictions with the existing experimental data. This 
model yields only a small and necessarily positive mag- 
netoconductance effect, whereas experimentally a large 
effect is found, either of positive or negative sign. It is 
found that whereas AEL/ EL is a first order effect, A I/I 
appears only in second order, in contradiction to the ex- 
perimental observation that they are of similar magni- 
tude. In addition, this model cannot account for the ex- 
perimental observation that AEL /EL measured at con- 
stant voltage greatly exceeds the same effect when mea- 
sured at constant current. In summary, we find that the 
excitonic pair mechanism model failed to explain OMAR, 
even only qualitatively. We trace the origin of the fail- 
ure of the model to the fact that it considers the spin- 
dynamics of polaron pairs, which, in first order, do not 
contribute to the current. 

By varying the injection efficiency for minority carriers 
in the devices, we show experimentally that AI /I is only 
weakly dependent on the EL quantum efficiency: AI / 1 oc 
t]yi with a ranging from 1/3 to 1/2. This dependence 
is unexpectedly weak if the effect were of an excitonic 
origin. This observation strengthens the conclusion from 
the modelling that OMAR is not due to an excitonic 
effect. 

Making the (presently unjustified) assumption that B 
acts directly on the carrier mobility, we found that the 
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qualitative features of both OMAR and the magnetic 
field effect on the electroluminescence naturally follow 
from this assumption. Therefore, we believe that a pair 
mechanism model acting on the current carriers directly 
is a promising direction for a future theory of OMAR. 
Spin-dependent bipolaron formation is a possible see- 



Next we calculate the expectation value of S z as a func- 
tion of time: 



S z (t)=<UU\S z \U1tl> (30) 
We obtain the following result (in units of h/2): 



VIII. APPENDIX: LORENTZIAN SHAPE OF 
MAGNETIC FIELD EFFECT TRACES DEDUCED 
FROM THE HYPERFINE HAMILTONIAN 

Here we will show that the empirical fitting formula 
AI/I oc B 2 /(B 2 + B 2 ) can be substantiated through a 
simple calculation to be presented in the following. The 
model considers the standard hyperfinc Hamiltonian, 



cos J uj 2 + a 2 t 



(31) 



Next we consider the case of a pair of spins each of 
which are subject to a separate Hamiltonian of form 
Equ. (|27|l . For simplicity we will treat the time-evolution 
of the initial state | ftT>i I 1U>2 only. Since the first 
spin is in an eigenstate it will not evolve with time, 
whereas the second spin's time evolution will be governed 
by Equ. I|29|) . We obtain therefore for the total spin (in 
units of h): 



H — u}qS z 



S I . 



(27) 



containing the electronic Zeeman energy and the hy- 
pcrfine interaction between a single electronic and nu- 
clear dipole. ojq — B , where g » 2 is the electronic 

g- factor and /is is the electronic Bohr magneton, S and 
/ are the electronic and nuclear spin (assumed to be 
1/2), respectively, and a is a measure (in units of fre- 
quency) of the hyperfine interaction strength. The z-axis 
is chosen to coincide with the direction of B. The Hamil- 
tonian will be written in matrix form where we use the 
following basis vectors: | 1*TT>, | JJ-T>) I 1U> an d | JJ4>- 
The boldface arrow denotes the z-component of the elec- 
tronic spin, whereas the second arrow denotes that of the 
nuclear spin. We obtain the following result: 



r i£o 

2 



H 



2 



-in 

2 







l£0 

2 



(28) 



It is evident from the form of H that | f|-|> and | JJ4> 
are eigenstates and therefore do not evolve with time 
other than through the trivial phase factor. However, 
| JJ.|> and | f|-j,> are mixed with each other through 
the off-diagonal matrix element. For simplicity we will 
now consider the time evolution of the | f|-J,> state only. 
It turns out that a calculation of the time evolution of 
the most general state vector and subsequent averaging 
over all electronic and nuclear spin orientations leads to 
similar results. To obtain the time evolution operator we 
perform a matrix exponentiation: 



Sl+ 2 (t) = 



1 1 / 

2 + 2 (pP 



PAP) 



(33) 



where pp and pap are the probability for finding the 
pair in a parallel or antiparallcl state, respectively. It is 
seen that Sl +2 (t) oscillates with time and that the peak- 



to-peak modulation depth is given by 



ui+a 2 



At large 



B, Sl +2 (t) remains close to 1 at all times although it 
performs a high-frequency (but small-amplitude) oscilla- 
tion. At small B, the frequency of the oscillation becomes 
smaller but its amplitude increases. The question arises 
whether the frequency or the amplitude of the oscilla- 
tion is the correct measure for the " spin- flip efficiency" . 
Because experiment shows (see section IIV A() that the 
oscillation frequency is much larger than the pair recom- 
bination rate, 7, it is the time average of S], +2 {t) that 
enters into the transition rate, specifically we may write: 



A7 

7 



PPlP + PAPlAP 

2 2(oj 2 + a 2 ) 
1 (B)- 1 (B = 0) 
7 (B = 0) 



1p + 



(34) 
"2T7ap(35) 



7-P — 1AP 



u>l + a 2 7p + ry A p 



(36) 



jp and 7ap are the recombination rates for parallel 
and antiparallel pairs, respectively. For a state initially 
in an antiparallel state 7p and 7^p have to be exchanged 
in Equ. ®. 

Finally we relate our results to the experimentally re- 
ported values, an — -£r~> f° r the hyperfine coupling 
strength: 



U = exp 



Ht 



(29) 



MFE oc 



B 2 



B 2 



(37) 



H 
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Schulten and coworkers [36( have shown that if the 
electron spin interacts with a large number of nuclear 
spins (as in the case of organic semiconductors) , then the 
individual an.i have to be added in a random-walk-likc 
manner. The final relation for Bq is therefore: 



1/2 



(38) 
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